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Abstract 

Results are presented of a fully ab-initio calculation of impact ionization rates 
in GaAs within the density functional theory framework, using a screened- 
exchange formalism and the highly precise all-electron full-potential linearized 
augmented plane wave (FLAPW) method. The calculated impact ionization 
rates show a marked orientation dependence in k space, indicating the strong 
restrictions imposed by the conservation of energy and momentum. This 
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anisotropy diminishes as the impacting electron energy increases. A Keldysh 

type fit performed on the energy-dependent rate shows a rather soft edge and 

a threshold energy greater than the direct band gap. The consistency with 

available Monte Carlo and empirical pseudopotential calculations shows the 

reliability of our approach and paves the way to ab-initio calculations of pair 

production rates in new and more complex materials. 
79.20.Ap,79.20.-m,71.15.Mb 
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Band-to-band impact ionization (I-I) is a carrier-carrier scattering process in which an 
energetic carrier creates an electron-hole pair through the excitation of a valence electron (e~) 
in the conduction bandSI This process is fundamental in small high-speed devices, both as 
a charge multiplication (e.g., avalanche photodiodes) and as a detrimental mechanism (e.g., 
field effect transistors). We present a fully first-principles approach based on screened- 
exchangei'S density functional theory (DFT)i that allows a full understanding of the basic 
mechanisms and physical quantities affecting the I-I process. Our final goal is to tune 
the rates according to technological requirements through band structure "engineering" for 
simple and complex materials, without need of ad-hoc parameters or pseudoatoms. 

All previous theoretical treatments have employed approximate band structures and ma- 
trix elements, based on k • p, Monte Carlo (MC)i~i or empirical pseudopotential (epp)HhI 
formalisms. To the best of our knowledge, the work herein reported is the first fully ab- 
initio calculation of I-I rates. We show results for the most studied direct gap semiconductor, 
GaAs. A variety of I-I results obtained through different approaches is available for GaAs, 
so that multiple comparisons can be made with results of the current method. 

We consider an e~-initiated I-I process (shown schematically in Fig. [I]). According to 
Fermi's Golden rule, the rate r(n 1 ,'k 1 ) at which the impacting e~ in a (ni,ki) state can 
produce I-I is obtained as: 
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where n 1 , n 2 , n 3 , n 4 are band indices (n^i = 1, 7 ) and k 1; k 2 , k 3 are k points in the full 
Brillouin zone (BZ). 5(E%? + — — -Eki+k 2 -k 3 ) shows the energy conservation. The 
antisymmetrized screened Coulomb matrix element is obtained by adding the probabilities 
in the singlet and triplet states \M\ 2 = \{\M D \ 2 + \M E \ 2 + \M D - M E \ 2 ) where M D and M E 
are direct and exchange (obtained from the direct term by exchanging final states) matrix 
elements. Each direct matrix element is expressed as: 
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^^(Ikx-ks + GuP + A 2 ) 
where the 5 function shows the momentum (k) conservation and Pn / ,k f ;n i) k i (G) is the Fourier 

transform of the overlap matrix of the wave functions, p n/j k f; ni,ki( r ) — ^n f k f ( r )^'ni,k i ( r )- The 
subscripts i and / denote initial and final states; e is the electron charge; Q is the unit cell 
volume; q = |ki — k 3 + Gu| is the momentum transfer and Go, Gu are reciprocal lattice 
vectors. 

The interaction between valence and conduction electrons is modeled by a Coulomb 
potential screened through a g-dependent static model dielectric function particu- 
larly accurate for semiconductors. The interaction between conduction electrons is mod- 
eled through a Debye potential, commonly used for the screening of impurity centers in 
semiconductors!. Both the temperature T and the carrier density in the conduction band 



n$ are taken into account through an inverse Debye screening length given by A = - / 47rn ° e 
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where Kb is the Boltzmann constant. Here, we used T = 300 K and no = 1-10 cm . To 
carefully evaluate the matrix elements, we used accurate wave functions and band structures. 
In particular, to overcome the well-known shortcomings of the local density approximation 
(LDA) to DFT when dealing with excited states!, we performed self-consistent screened ex- 
change (sX-LDA)!&0'0 calculations, as implemented within the highly precise full-potential 
linearized augmented plane wave (FLAPW)0 method. During the initial LDA (and the fol- 
lowing sX-LDA) self-consistent iterations, we used a cut-off equal to 9 Ry (6 Ry) in the 
wave vectors and I <8 (I <4) inside the muffin-tin spheres chosen as R^t = ^mt = 2-3 
a.u.. The summation over the irreducible BZ was done using 4 special k points@. 

As shown by Seidl et a/.!, the sX-LDA approach can be recast within a generalized DFT 
formalism in which the inclusion of a non-local sX functional highly improves the description 
of the conduction band states compared with a bare LDA approach. This improvement is 
essential in the present context, since the transitions considered always involve conduction 



states. The many-fold integration over the full BZ was carried out using the technique 
proposed by Sano and Yoshiii, based on a regular grid of k-points (with an interval length 



approximated by a "top-hat" function, i.e. as a rectangle Ae wide and 1/A^ high (we 
used A^ = 0.2 eV0). Spin-orbit coupling and Umklapp processes (Go 7^ 0) have been 
neglected; their inclusion in the formalism is the subject of ongoing work. We also have not 
considered: i) phonon-assisted transitions, that would relax the k conservation requirements 
among the four involved electronic states (i.e. k conservation would be satisfied through 
phonon participation) and it) the influence of the electric field on the collision term, i. e. the 
"intra-collisional field effect "0. 

The results for I-I rates initiated by electrons in the second lowest conduction band with 
wave vectors along [001] and [111] are shown in Fig.[|. The insets show the band structure 
along the same symmetry lines as in the main panels. Due to energy and k conservation 
constraints, there is a marked k-space anisotropy in the I-I rate, that is orders of magnitude 
higher along [111] (r to L) than along [001] (r to X). It is noted that the second lowest 
conduction band in GaAs shows a decrease in energy from T toward X, but increases with 
increasing |k| in other directions, such as [111]. We therefore expect, and find (see Fig.^](a)), 
a "wave-vector anti-threshold" along the [001] axis at which I-I no longer becomes possible. 
Our calculated points are compared with those of other calculations^ obtained with EPP 
that investigated the effect on the predicted rate assuming two different band structures0'0. 
This comparison makes clear the importance of employing an accurate band structure (such 
as sX-LDA FLAPW), to obtain reliable rates. The degree of agreement between these 
calculations and ours is reasonable, given that the band structure and some of the numerical 
approximations^ made in the evaluation of matrix elements differ. 

The calculated I-I rates for GaAs are shown in Fig.|3| (a) for processes initiated by elec- 
trons in the three lowest conduction bands vs. their impacting energy. The scattered points 
in the low-energy region reflect the strong anisotropy already noted: carriers with the same 
impacting energy but different wave vectors can have widely varying rates. However, this 




As in previous investigations^'!, the energy 5 function i 
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anisotropy diminishes at higher energies, due to the greater ease with which k and energy 
conservation restrictions can be satisfied. Moreover, it is of interest to show an "isotropic" 
I-I rate, that depends only on the impacting electron energy E (solid line in Fig. [3] (a)): 

The physical reason to obtain an isotropic energy dependent rather than a wave-vector 
dependent rate is that in most common technological devices in which high electric fields 
usually are present, carriers are scattered by phonons, so as to reach similar energies, but 
largely different wavevectors. 

Further, our independent rate has been fitted using a "Keldysh-type" formula (dashed 
line in Fig. |3] (a)), R(E) = P[E — E t h] a where E t h is the "isotropic" threshold energy, 
i.e. the minimum energy at which the carrier is able to excite an e~ in the valence band 
and, therefore, to initiate I-I. Here, P, a and E t h have been treated as fitting parameters. 
An optimized linear regression procedure yields a fitted value E{^ ~ 1.8 eV, whereas our 
simulations indicate no ionization events for impacting energies lower than E^ lc ~ 1.86 eV 
(considered as our "real" value for Eth)- The excellent agreement between E(^ and E^ lc 
provides confidence in the numerical fit. As a result, the threshold is slightly higher than the 
energy gap and in good agreement with the relation E th = 3/2E gap obtained from a parabolic 
band structure with constant effective massesEI. However, we believe this agreement to be 
fortuitous, since we have shown previously the importance of a careful treatment of band 
anisotropy. Moreover, it is noted that the value of E t h is exact within the uncertainty given 
by From the fit, we obtained P = 3.5xlO w s~ 1 eV~ a and a ~ 5.8; the high value of a 
reveals the "soft" character of the GaAs threshold. 

In Table | and Fig. |3] (b) we compare the result of our fit with other Keldysh-type fits 
available from M(J or EPP00 methods. Our ab-initio results are in overall good agreement 
- especially at low energies - with two of the previous workffl. While the threshold energies 
are within a very similar range, there are some differences in the P and a values; these can be 
ascribed mainly to details of band structures and numerical methods employed for evaluating 



the rates (e.g. many-fold integration schemes). 

While it would be important to compare our results with experiment, the I-I rate is not 
directly comparable. Rather, the quantity commonly measured is the electron ionization 
coefficient a(F) (in cm -1 ) as a function of the applied field F. This coefficient is related to 
the inverse of the mean distance traveled by carrier prior to I-I; contrarily, our calculated 
rates (when summed with the phonon rates and normalized) give the probability (in s _1 ) 
of an I-I event as a function of the carrier's energy. Therefore, the calculated I-I rates and 
the measured I-I coefficient are not easily related. For the experiments, we refer to Bulman 
et a/.S, who measured a(F) in (100) GaAs in a large number of different p + n structures 
using avalanche noise and photocurrent multiplication: their electron ionization threshold is 
1.7 eV. Therefore, the agreement with our value of ~ 1.8 eV supports the reliability of our 
procedure. 

Finally, we offer some information about the distribution of final states, that is of great 
interest for transport simulations. In Fig. ^ (a), we plot the average final electron energy 
vs impacting electron energy along with the best linear fit. Similar to the rates themselves, 
the scattering of points off the straight line is evident in the low energy region. The average 
final electron energy is not necessarily equal, even for primary impacting electrons having 
the same initial energy, because of the strong restrictions imposed by k conservation. On the 
other hand, as the impacting energy increases the linear fit improves significantly. Moreover, 
we plot in Fig.§ (b) the percentage of transitions that involve one or both final states in the 
symmetry point valleys (L, T, X) vs the impacting electron energies. As expected, at low 
energy, most (even 100 % at some energies) of the transitions involve one or even both final 
states in the T valley. However, as the impacting electron energy increases, the other valleys 
(especially the L valley) become accessible and the final states are generally more spread 
over the BZ. 
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TABLES 

TABLE I. Fitting parameters for the rates shown in Fig. 3 (b), using the Keldysh fit formula 
R{E) = P[E-E th f. 
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FIGURES 

FIG. 1. Electron- initiated I-I process. The initial electrons in the conduction and valence bands 
are in states 1 and 2, respectively; after the transition, the final electrons in the conduction bands 
are in states 3 and 4. 

FIG. 2. Calculated ionization rates (in s _1 ) for impacting e~ in the second conduction band 
and wave vector along ((a) T — X and (b) T — L) vs normalized wave vector (filled diamonds). 
The open symbols show previous results from Ref. [H]. In each panel, the inset shows the band 
structure along the corresponding symmetry lines (the zero of energy is set to the conduction band 
minimum (CBM)). 

FIG. 3. (a) Ionization rates (in s _1 ) for initial e~ in the three lowest conduction bands vs their 
energies (the zero of energy is set to the CBM). Filled circles: wave vector dependent ionization rate 
.R(k); solid line: energy dependent ionization rate R(E); dashed line: Keldysh fit. (b) Keldysh-type 



fits taken from this work (dashed bold line), Ref. |10| (thin solid line), Ref. 11 (dot-dashed line) and 
Ref. @ (dotted line). 

FIG. 4. Panel (a): Average final e~ energy (filled circles) vs impacting e~ energy (in eV) 
referred to the conduction band minimum (zero of energy). The solid line shows the best linear fit. 
Panel (b): Percentage of transitions that involve final states at the symmetry point (r (triangles), 
L (diamonds) and X (circles)) valleys vs impacting e~ energies (in eV). 
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